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1. 


Introduction . 


In  the  single-server  queueing  iTiOdel  with  uniformly  bounded 
virtual  waiting  times  studied  here, the  server  never  permits  the 
amount  of  work  still  to  be  handled  by  him  at  any  moment  to  exceed 
a  constant  D.  This  means  that  the  total  time  spent  by  each 
customer  in  the  system  is  limited  by  D.  If  this  time  exceeds  D, 
then  the  customer  has  to  cut  short  his  assigned  service  time. 

This  model  is  also  the  model  for  the  dam  with  finite  capacity, 
instantaneous  water  supply  and  constant  release  rule.  The  content 
of  the  dam  is  equivalent  to  the  virtual  waiting  time  of  t.he  queueing 
system.  If  at  some  input  moment,  the  supply  is  so  large  that  the 
content  exceeds  the  capacity,  then  overflow  occurs. 

The  Ej/D/1  model  has  been  investigated  by  Ghosal  [9]  and 
Prabhu  [15];  the  M/G/1  model  by  Kovalenko  [11],  Takacs  [23]  and 
Cohen  [5].  Cohen  [4]  also  obtained  several  results  for  the  model 
in  which  the  inter-arrival  times  and  service  times  have  rational 
Laplace-Stiel jes  transforms.  Daley  [7,3]  studied  the  distribution 
of  the  stationary  waiting  time  of  the  GI/G/1  model.  (See  also  [1] 
and  [2]).  In  this  paper,  we  shall  obtain  a  fairly  complete 
solution  for  the  GI/G/1  model.  Explicit  results  will  be  obtained 
for  the  M/G/1  model. 

The  key  to  our  analysis  of  the  system  is  that  many  of  its 
processes  are  regenerative ;  that  is,  they  restart  probabilistically 
whenever  a  customer  initiates  a  busy  period.  Regenerative  processes 
in  this  sense  were  introduced  by  Smith  [l7,18]  and  have  been  used 
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by  many  authors  to  study  the  stationary  behaviour  of  many 
queueing  systems.  (  See,  for  example,  [3],  [5],  [$ ]  ,  [12],  ll3  '  , 
[16]  [19]  )  .  By  using  the  regenerative  properties  along  with 

analytical  methods  in  this  paper,  we  shall  show  that  not  only 
the  stationary  behaviour  of  the  system  can  be  studied  but  its 
transient  characteristics  can  also  be  obtai.ned.  These  methods 
also  give  us  insight  into  the  probabilistic  structure  of  the 
system.  In  Section  3,  the  mathematical  description  of  the 
transient  behaviour  of  the  system  is  obtained  from  its  behaviour 
within  a  busy  cycle.  In  Section  4,  the  mathematical  description 
of  its  stationary  behaviour  is  also  obtained  from  its  behaviour 
within  a  busy  cycle.  The  behaviour  of  the  system  within  a  busy 
cycle  and  t.he  stochastic  laws  for  the  busy  cycle  are  then 
studied  in  Section  5.  In  Section  6,  due  to  the  special  structure 
of  the  M/G/1  queue,  we  shall  obtain  explicit  results  for  the 
behaviour  within  a  busy  cycle  and  the  stochastic  laws  for  the 
busy  cycle  of  this  system. 
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2 .  The  formal  model  and  notation. 

We  are  given 

(D.l)  a  real,  positive  number  D,  D  <  *; 

(D.2)  an  integer-valued,  non-negative  random  variable  m^ , 

(0.3)  a  real,  non-negative  random  variable  w^,  P’"wq<  D}  =  1; 

(D.4)  Two  independent  sequences  of  independent  and  identically 
distributed,  real,  positive  random  variables  {tj^,k>  1}  and 
{sj^,  k>  1}.  We  assume  that  each  of  and  t^  has  a  finite  first 
moment. 


m^  is  the  number  of  customers  in  the  system  at  time  t=0- 

and  Wq  is  the  virtual  waiting  time  at  time  t=0-.  Let  customers 

(m„+l) ,  (m.+2) ,  ...,  k  arrive  at  the  time  epochs  t  ,  ,, 

-u  -u  ^ 

irngt  2 . Ik'  •••  0  =  Im^+l  "  ImQ+2  "•••"  ^k  " 

Let  “  Ik  ~  -k  ' 

Let  the  assigned  service  time  of  the  k’"  customer  be  Sj^. 
We  assume  that  no  customer  can  stay  in  the  system  longer  than  the 
time  interval  of  length  D.  If  the  total  actual  waiting  time  and 


service  time  of  a  customer  exceeds  D,  then  he  shall  cut  short  his 


assigned  service  time. 


We  write 

(D.5)  (z)  ”  E  {exp  (-zSj^)  }  fgj-  Re(z)  >  Or 
(D.6)  f2(z)  =  Etexp(-ztj^)  }  for  Re(z)  >  0. 


The  customers  are  served  in  order  of  their  arrivals  and  there  is  no 
limit  on  the  size  of  the  waiting  room. 
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(  D  .  7  ) 
(D.  8j 
(D.9) 

(D.IO) 
(D.lx) 
(D. 12) 
(D.13) 
(D.14) 
(D.15) 
(D.16) 
(D.17) 


(0.18) 

(D.19) 

(D.20) 

(D.21) 


(D.22) 


Besides  t,  ,  we  want  to  study  the  followinc  random  variables: 

-K 

w,^  =  the  actual  waitin'  time  of  the  k  ‘  customer  Oom^^)  ; 

w  ,  =  limit  in  distribution  of  w  when  k^>-',  if  this  exists; 

—  X*  K 

i.,  =  the  lost  service  time  of  the  k^^  customer  (k>mp^)  , 
i.e.  =  niax  (Wj^  +  Sj^,D)  -  D 
.1  =  limit  in  distribution  of  i  when  k-^“,  if  this  exists; 

=  the  duration  of  the  initial  busy  period; 
g^  =  the  duration  of  the  busy  period,  v>  2; 

i^  =  the  duration  of  the  first  idle  period; 
i^  =  the  duration  of  the  idle  period,  v>  2; 

Cj^  =  g^  +  ij^  =  the  duration  of  the  initial  busy  cycle; 
c  =  p  -t-  i  =  the  duration  of  the  busy  cycle,  v>  2; 

=  the  number  of  customers  served  during  the  initial  busy 
period,  including  the  customers  in  the  system  at 
time  t=0-; 

t  h 

n^  =  the  number  of  customers  served  during  the  v^  busy  period, 
v>  2; 

v{t)  =  the  virtual  waiting  time  at  time  t,  t>  0; 
v(”)  =  limit  in  distribution  of  v(t)  when  t-^°°,  if  this  exists; 
)^{t)  =  the  total  number  of  customers  arriving  during  the  time 
interval  [0,t],  including  the  m^  customers  in  the 
system  at  time  t=0-; 

a(t)  =  the  time  difference  between  t  and  the  time  of  the  first 
arrival  during  the  interval  (t,“) . 


The  results  will  be  expressed  in  the  following  forms: 
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(D.23)  W^_„(x,^,z)  =  E-:exp(-5T^-zw^)imQ=m,WQ=w} 


k  ‘■rik'  ‘-0  “'-0 


for  0<!x|<l,  Re(^)>0,  'z<=°,  m>0 ,  w^O  ; 


=  Ik=m+l^^E{exp(-Cl3^-2ik)  lmQ=m,WQ=w} 

for  0<|x|<l,  Re(C)>0,  Re(z)^0,  m>0 ,  w>0; 

*CO 

(D.25)  V  (x,^,z,s)  =  exp(-Ct)E{x-^^Uxp{-zv  (t) -sa  (t)  )  I  rn  rm,w  -w)dt 
in ,  w  Jq  ^  — <j- 

for  0<lx|<l,  Re{C)>0,  | z j <«, Re (s) >0 ,  m>0,  w>0; 

-1 

{D.26)  ^(x,C,-z)  =  E^x  exp(-5Cj^+z^^)  |mQ=m,WQ=w} 

for  0<|x|<l,  Re(5)5'0,  Re{z)^0,  m^O ,  w^O; 


{D.27) 

II 

N 

8 

£{exp(-zY^)  } 

for 

1 z 1 <“  ; 

(D.28) 

Loo(z)  = 

Ei'exp  (-2^^)  } 

for 

Re(z)>0  ; 

{D.29) 

II 

N 

8 

> 

E{exp (-ZV (») } 

for 

1  z  1  <» 

We  shall  need  the  following  intermediate  Laplace-Stiel jes 
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It  is  important  to  point  out  here  that  there  are  two  types 

of  busy  periods  which  Cohen  [  4,p.284]  calls  strong  and  weak  busy 

periods.  While  two  consecutive  strong  busy  periods  are  separated 

with  probability  one  by  an  idle  period  of  non-zero  duration;  a  weak 

busy  period  may  be  followed  by  an  idle  period  of  zero  duration. 

t  hi 

In  other  words,  if  the  (k+1)  customer  arrives  at  t.he  instant 
the  k^^  customer  departs,  the  strong  busy  period  continues  while 
the  weak  busy  period  terminates  and  a  new  weak  busy  period  starts. 
We  consider  both  types  of  busy  periods  in  this  paper.  If  the 
result  is  applicable  for  both,  no  notational  distinctions  are 
made.  If  a  result  is  applicable  to  the  strong  busy  period  only, 
then  a  superscript  "s"  is  added  to  the  notation.  If  a  result  is 
applicable  to  the  weak  busy  period  only,  then  the  superscript  ”w" 
is  added. 

Remarks ; 

(R.l)  From  (D.?0)  and  (D.31),  we  have 

(2.1)  E {n, I m-=m, w„=w}  =  W  {l,0,0)+m  =  L_  (l,0,0)+m  ; 

— 1  — 0  —0  m ,  w  ni ,  w 

(2.2)  E^Il2^  “  ^  ^0,0 


(1,0,0) 


3. 


Regenerative  results  for  the  transient  behaviour. 


In  this  section,  we  shall  show  that  the  study  of  the  actual 
waiting  time  of  each  customer,  the  lost  service  time  of  each 
customer  and  the  virtual  waiting  time  at  each  epoch  can  be 
reduced  to  the  study  of  the  lost  servive  times  within  one  busy 
cycle  alone;  that  is,  for  0<|x|<l,  Re (5)-  Re(s)sO,  m>0,  wsO, 
W^^(x,^,z)  (  I  z  !  <«)  ,L^  ^(x,C  ,z)  (Re(z)iC)  and  ^  (x ,  C  ,  z  ,  s)  (;z'<=°) 

can  be  obtained  from  L  (x,f,2)  (Re{z;  ^0)  and  C  (x,^,0). 

m ,  w  m ,  w 

The  arguments  are  based  mainly  on  the  regenerative  property 
that  the  continuations  of  many  processes  in  this  system  beyond 
the  end  of  a  busy  cycle  are  the  probabilistic  replicas  of  these 
processes  commencing  at  the  beginning  of  that  busy  cycle. 


THEOREM 

_1  For  0< 

i  X 1 < 1 ,  Re ( C)  >0 , 

1  z  1  <'» 

(3.1) 

=  Wq^q(X,C,2)  / 

[1-C 

(3.2) 

W  (x,c,z) 

m,w 

=  W  (x,C,z)  + 

m ,  w 

C  ( 
m,w 

PROOF : 

We  have  , 

from  (D.23)  and 

(D.30) 

(3.3) 

W  (x,C,z) 

m,  w 

=  W  (x,C,z) 

m,w 

0,0 


0,0 


+  E- 


{S.l*  ^  +k-"!in,+k> 


-1 


mQ=m , Wq=w 


Now  since  the  queueing  process  starts  again  probabilistically 
when  the  (n^+1) customer  initiates  the  second  busy  period  at 
time  T  =  c,  ,  c,  and  n,  are  independent  of  w„  for  all 

— iij^+l  —1  —1  —1 

k.^1.  Also,  for  all  k^il,  w^  are  independent  of  and 


Thus  we  can  write 


(3.4) 


'•J  .,(x,^,z) 

m ,  A 


t“k=l 


X  exD  {  -  ^T,  -zw 


T, -zw  im-  =  0,'A-  =0^ 

— k  — n, +ki— 0  —0  J 


Furthermore,  if  ^=0  and  ^=0,then  w  ,,  will  have  the 
— u  — u 

same  distribution  as  w^^.  Thus  (3.2)  is  proved.  As  (3.2) 
is  also  applicable  when  51^=0,  iiQ=0  /  (3.1)  follows. I 


THEOREM  2 ;  For  0<|X|'<1,  Re(:^)>0,  Re(z)>0,  m>  0 ,  w>0, 

(3.5)  LQ^g(x,C,z)  =  L.^q(x,C,z)  /  [i-CQ^Q(x,;,0)]] 


(3.6) 


,  (x  ,  ^ , 0 ) Lq  q  ( ^ 


PROOF ;  The  proof  is  omitted  beca\ise  it  isi  analogous  to  that  of 
Theorem  1.1 


Theorems  1  and  2  relate  W  (x,C/Z)  and  L  (x,C/Z)  to 

m, w  m ,w 

A  A 

W  (x,^,z),  L  (x,CfZ)  and  C  (x,^,0).  We  shall  now  show 
m  I  w  rn  /  w  m  ,w 

that  ^(x,C,z,s)  can  be  obtained  from  ^(x,C,z)  and 

^m  This  is  an  important  relation  which  is  of  interest 

by  itself  because  it  enables  us  to  find  the  mathematical  description 

of  the  behaviour  of  a  queue  in  continuous  time  if  we  know  its 

behaviour  at  a  certain  set  of  discrete-time  epochs. 

LEMMA  1;  For  0<(xf<l,  Re(C)>0,  | z ( <«  ,  m>0,  w>0, 

(3.7)  Vq^q(x,5,z,s)  =  Vq^q(x,C,z,s)  /  [l-CQ^Q(x,C,oj] 


PROOF ;  The  proof  of  this  lemma  is  omitted  because  it  is  analogous 


to  that  '-f  '^leorc-.'ii  1.  Here,  we  u.^e  th 


.  »  '  -I  -  , 

‘-1,  \  / 


:hien  icr  all  t/.O,  ‘  - 


same  ais 


distr  ib'at  ic;-. :s  ■  j  ( • 


:  I  an  :  i 


LEMMA  2;  For  0'!x'rl,  Re  {  H  >.  Re  (  s)  > 


(3.9)  (2  +  s-=:)  l,z,s) 


,  n  ■  U  ,  v;  •  0  , 


-x™exr(-zw)  -  zx^'exp  (  (s- 1)  w)  /  ( s- 1 ) 

r  n  - 

+  I  l-x3’(z)Q(s) 

r  1  ^ 

t  z  i-xt (^-s) n(s)jw^^,^(x, ',t:-s)/x(s--;) 

+  exp(-zD)  P.(s)  [Lj^^^(x,-:,z)-Lj^^^_^(x,';,h) 


PROOF : 


+  zexp  ( (s-C)  D)  r.’ (s)  1 ,  l-s) (x,  ,  0)j  /(s--,) 

+  (  z  +  s - 1 ,  ( X ,  C , 0 ) / ( s - 1 ) 

n ,  w 

For  U<|xUl,  RG(^)^0,  Re(s)3.0,  |z!<«  ,  -1 )  ,  w>0, 


(3.10)  V'  (x,t:,z,s)  = 

rn  f  w 


-jt-l  exp(-Ct-2{T  +min(w^+Sj^.D) -t) )  dtlm^  m,<^ 

-\^k=iTi+l  JX)^  ^ 

e{x"M^  exp(-5t-x(l„  +mln(w  ,D) -t) -s 

I  In,  ^  ^ 


+  exp(-5t-s(T^  +l-t))dt|niQ=in,^-wj 

I  El  —1 

AS  =  ik  '  -k+i 


n  =  T  +  inin(w  +s  ,D)  ,  (3.10)  becomes 
^1  -nj^  -El  -Hi 


for  k<nj^  and 
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1  nJr-r^  kT  ,  . 

Z+S-?  E{^k=m-1  ^  L^^P(-^ik+l-^^k+l’ 

-  exp(-lT^-2(:nin(Wj^+Sj^,D)  )-stj^)]  jmQ=m,w^=w  | 

1  )'  -1  r 

^  iTiTr  e|  X  ^exp{-S£^-s(c^-D^)) 

-  exp(-CT^^-z(min(w^^^s^^,D))-st^^)]|n^=^  | 

s-C  ■'  ~  j^xp(-'c^)  -  (c^-£^)  )~j  |ii[^=in,^z=v^  I 

Now  observe  that  for  k^n^ ,  =  max  (Wj^+s^,  D) -D  and  hence 

(3.12)  exp(-CXj^-z  (min  (Wj^+Sj^,D)  )-stj^) 

=  ®XP  (W)^-^Sj^) -stj^)  +  exp  (-C2j^“2D-stj^) 

-  exp(-CT^-z  (max(Wj^+Sj^,D) ) -stj,) 

=  exp  ( -  z  s )  exp  ( -  s tj^)  exp  ( -  C  T  zw^ ) 

-  exp(-zD)exp(-stj^)  j^xp(-5T^-zi^) -exp(-?T_j^)~| 

Thus  (3.11)  can  now  be  written  as 

.A 

=-x^exp  (-zw)  /  (z+s-^ 

+  [l-x'i'(z)S^(s)JWj^^^(x,5,z)  /  x(z+s-5) 

+  exp(-zD)f2(s)[L^^^(x,C,z)  -  L^,^(x,C,oJ  /  (z+s-C) 

/  (z+s-C) 


+ 


/  (s-C) 


It  is  easy  to  prove  that  ,3.^3'  is  also  applioable 

when  n,  =  ni+i .  xow  since  V  (x,3,z,s)  ,?  (x,-:,3)  and 

""-L  IT.  /  w  m ,  V 

(x,4/0)  are  analytic  -^cr  0<  x  s;l.  Red)  ^  Re(s)^C,  z,<^f 
letting  z='-s  in  (3.1^)  yields 

(3.14)  ^(x,-;,s)  =  x^exp(  (s-£)  w) 

-  i  l-xf  (C-s)  r.{s)  ,W^^^(xdd-s)/x 

If  we  substitute  this  equation  back  into  (3.13)  ,  we  obtain  (3.9  ) 

THEOREM  3  :  For  0<;x|<l,  Re  (  ^)  _^e  ( s)  ^0  ,  \z]<x,  m^Q  ,  w^O  , 

(3.15)  (z+s-OV  (x,l,z,s)  =  -x”'exo{-zw)  -  zx^exp  (  (s-dw)/(s 

m ,  w 

+  j^-x4' (z)n(s)  jw^  ^(x,C»z>/x 

+  z[T-x'V  (C-s)  Cl(s)  tw  (x,;  d-s) /x(s-t) 

+  exp  (-20)  (s)  ^(xddTl 

+  zexr)(  (s-' )  D)n  (s)  fL_  (x,'d-s)-L„  ix  ,  ^  ,  oT\  /  (s-’^) 

PROOF ;  The  proof  is  straightforward  from  Theorems  1,2  and 
Lemma  s  1,2.1 

It  remains  to  show  that  W  (x,5,z)  can  be  obtained  from 

m,  w 

L  (x,l,z)  and  C  (x,C,-z). 
m ,  w  m ,  w 

THEOREM  4 ;  For  0<|xi^l,  Re  (?) ^Re (z) >0 ,m^0 ,  w^O , 

(3.16)  [i-x'(' (z)  n  (^-2)  ]  W_  „(xd  ,  z) 

in  /  w 

=  x”'"*’^exp  (-zw)  -  xC  (x,5,-2) 

■■  m ,  w 

-  xexp(-2D)n('-2)  [Lj^  .^(x,f  ,2) -L^^^(x,d0)  - 


Remarks : 

(R.2)  (3.1)  and  (3.7)  are  the  generalizations  of  (3.3!  and 

(4.6)  in  [5/Pps.  6,131  respectively. 

(R.3)  When  ,  then  Theorem  3  becomes  Theorem  2  in  121]  . 

While  Takacs  derived  the  'latter  directly,  the  former  is  obtaine 
via  Lemma  2,  which  will  also  be  useful  in  the  derivation  of 
Theorem  6  later. 

(R.4)  If  we  let  x^l  ,  z-^0,  s->-0  in  (3.14)  and  (3.16)  and  then 
use  1' Hospital's  Rule  to  obtain  the  limit  when  i-'O,  we  shall 
obtain  the  following  Generalized  Wald's  Lemma: 


(3.17)  £{£!> 

(3.18)  e:£^} 


E{Wq}  +  "Ej lk=m^i:^k;  ' 

[E{n_j_}-E{mQ})E{t^}  • 


has  beer. 


4  . 


Regenerative  results  for  the  stationary  behavioar. 


For  the  queueing  system  studied  in  rhis  paper,  ir 
proved  that,  unless  ~  ~  consta.nt  <  D,  each  of  the  processes 

{Wj^,k>mo},  {v(t),t>0;  has  a  unique  stationary 

distribution  which  is  independent  of  the  initial  conditions 
(See  [7] ),  In  this  section,  we  shall  show  that  the  study  of  the 
stationary  behaviour  of  the  system  can  also  be  reduced  to  the 
study  of  the  lost  service  times  within  one  busy  cycle  alone;  that 
is,  the  expressions  for  W^(z)  (!z'<^),  L^(z)  (Re{z)^0)  and 

V^(z)  (|zl<‘p)  can  be  obtained  from  Lq  q{1,0,z).  Here,  we  shall 

use  a  general  theorem  in  the  literature  stating  that  the  stationary 
distribution  of  a  regenerative  process,  if  it  exists,  is  the 
'time  average'  or  'customer  average'  of  the  process  over  a  regene¬ 
rative  cycle.  (See  [  5J) 

This  allows  us  to  state  the  following  theorem  without  proof : 


THEOREM  5: 


(4 

.1) 

W  ( z ) 

00 

Wg,o(l,0,z)/E{n2} 

for  i  z  1  <=o 

(4 

.2) 

II 

N 

8 

hg 

for  Re(z)^0 

Together  with  (3.16),  (2.2)-,  the  assertion  for  W  (z)  and  L  (z) 
is  now  true.  The  next  theorem  will  allow  us  to  find  the  distribution 
function  of  the  stationary  virtual  waiting  time  in  terms  of  the 
distrioution  functions  of  the  stationary  actual  waiting  time  and 
the  stationary  lost  service  time. 
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THEOREM  6 ;  For  z 

(4.3)  V^(z)  =  1  -  E '-£2 {1-r' (z)  (z) /zE'.'t^; 

+  exp (-zD)  [L^ (z) -l]/zE{t^ ; 

PROOF:  This  is  because  V  (z)  =  V-  ^  ( 1 , 0 , z , 0 ) /£ { c_ }  .  Upon 
applying  I'Hospital  rule  to  (3.13),  we  prove  the  theorem.  £J 

Remarks : 

(R.5)  When  ,  then  (4.3)  beccfraes  a  well-known  result  due  to 
Tak/cs  [20.]  for  the  classical  GI/G/1  queue. 


Ths  stcchastic  laws  for  the  busv  cvcles. 


Let 

(D.34)  u  =  the  set  of  all  those  functions  of  z  which  are  analvtic 

I  <2 

in  the  domain  Re(2)>0  and  continuous,  free  from  zeros, 
'uniformly  bounded  in  Re(z}>0; 

(D.35)  =s  the  set  of  all  those  fvinctions  of  z  which  are  analvtic 

"Z 

in  the  domain  Re(z)<0  and  continuous,  free  from  zero^ , 
uniformly  bounded  in  Re(z)<0; 

(D.35)  D  =  the  set  of  all  those  functions  <fr(z)  which  are  defined 
z 

for  Re{2)=0  on  the  complex  plane  and  can  be  represented 
in  the  form 


(5.1) 


^(z)  =  E{£exp(-zn)}  , 


where  £  is  a  complex  (or  real)  random  variable  with 
and  ■  £  is  a  real  random  variable. 


Let  us  define  the  following  transformations  on  r^: 


(D.37) 

(z)  } 

= 

EC^exp  (-zn"^)  } 

(D.38) 

T*{'f  (2)  } 

3 

*^(2)  -  T^{^(Z)}  ; 

(D.39) 

U„{'P  (2)  } 

=» 

E{£6  (n2  0)exp  (-zji)  } 

(D. 40) 

(z) } 

= 

^  (z)  -  U^(^lz)}=  eU- (£<0)exp(~zr) }  ; 

(D.41) 

V^(f(z)} 

3 

E{£3 (n>0)exp (-zn) } 

(D.42) 

3 

^(z)  -  =  EilUr-O)  exp(-z-)  }  , 

where  n"'’ 

=  max ( 0 , n ) 

and 

6  (A)  is  the  indicator  fxinction  of  any 

i!* 

event  A; 

that  is ,  5 (A) 

=1  i 

r  A  occurs  and  6(A)=0  if  A  does  not  occur. 

Clearly,  = 


inc  T*‘^'^(  =  '  ■'  U*';Mz)}  and  belon^j  tc 

2  2  •" 


.  Aisc,  It 


is  easy  to  shew  that 

(5.2)  U_{?(z)}  =  T,{5(z)}  +  lim  T*{J)(z)’ 

2  2^^®  '  2  ' 

(5.3)  U*{5(z)}  =  -  lin  7*{>[z)} 

£t  2  2^*®®  *  2  '  ' 

k 

(5.4)  v,{^(z)}  =  ' 

(5.5)  v*{^’(z)}  =  T*^5(z)}  +  T2^ (2)  }  ; 

(5.6)  T20(z)}  =  U2^*(2)}  +  11^2-0 

(5.7)  T2f'*’(z)}  =  V2^'^(2)}  +  lira^_^Q  yH'^iz)) 

This  means  that  the  closed  form  expressions  for  these 
transformations  can  be  obtained  if  that  for  known. 


The  following  lemma,  which  is  due  to 

enable  us  to  obtain  T-t^’fz)}  explicitly: 

z 


Taj<acs  [22  J  ,  will 


LEMMA  3:  If  4>(z)«R_  ,  then  for  Re(z)>0,  we  have 

”  z 

(5.8)  =  HHO)  +  lim^^Q  ^  jlll^s  , 

where  the  path  of  integration  (oO)  consists  of  the  imaginary 
axis  from  z=-i«  to  2=-i£  and  again  from  z*ie  to  z*i*. 

PROOF:  See  Theorem  2  in  [22  3.0 


In  this  paper,  we  have  shown  in  Sections  3  and  4  that  both 


the  transient  and  stationary  behaviours  of  the  system  can  be  studied 
in  terra  of  the  lost  service  times  within  one  busy  cycle  alone. 


In  this  sectitn ,  analytic  r.ethcds  are  given  fcr  finding  the 

integral  ec’iations  that  would  theoretically  allow  us  to  obtain 

results  for  the  lost  service  tirtes  within  a  busv  cvcle  1  z: 

"  *  .t ,  w 

and  the  stochastic  laws  of  the  busv  cvcle  C  -z)  si.T.ultaneous 

m ,  w  ’ 

These  equations  will  be  expressed  in  terms  of  the  transformations 
defined  in  (D . 37) - (D . 4 2) . 

Basically,  this  method  simply  involves  the  re-arrangem.ent  of 
(3.16)  into  identities  whose  left  hand  sides  belong  to  and 

right  hand  sides  belong  to  .  By  Liouville's  Theorem,  they  are 
functions  independent  of  z.  The  integral  equations  will  be  obvious 
when  these  functions  are  known. 

First,  for  the  sake  of  simplicity,  let  us  write 


for  0<lxl<l,  Re(C)>0,  Re(z)^0,  m>0 ,  w^O. 


From  (3.16),  we  have 


(5.9)  L  ,  (x,;,0)  =  „(x,?,oy  =  Ix^  -xC  „(x,5,0)  ]/[l-xri(C)  ]  . 

m,w  ra,w  m,w 


This  means  that  L  (x,£,z)  will  be  known  if  Q  (x,£,z)  and 

m ,  w  ^  m ,  w  ^ 

C  (x,5,-z)  are  known, 
m,  w 

For  0<|x|<l,  O^Re (z)^^^ (£) ,  we  now  assert  that 
[!-■!' (z)  P.  (C-z)  ]  can  be  factorized  into  the  form 

(5.10)  [l-x^*  (z)f2  (C-z)  ]  =  g''’(x,C#z)/g”(x,Cf  z)  , 

where  g'*’ (x , '  ,  z)  €  g~(x,C,z)  ^  n^. 

Such  factorization  always  exists  as  we  can  write 

(5.11)  g'*‘{x,5/Z)  =  expjT^dnd-x'i' (z)n  (--z)  )  }| 

-  expj-Tjxln  (1-x'y  (z)n  (i-z)  )  }l 


(5.12) 


g  (x,5,z) 


Ir.  fac”,  ::  (x,^,Z/  and  5  are  dererr.ir.ed  up  zo  a  rr.ulriplucat 

function  of  x  and  For  if  we  also  have  '  1-x  r  ( 2  }  i ; -z )  '  = 

h~  (x  ,  ^  ,  2; /"r.”  ;x  , ,  2  )  where  h  (x,^,2)  =  ;x,;,2)  '  sj  ^ 

Chen  bp-  Liouville's  Theorem  g""  (x,  ^  ,  2 ) /h  '  (x ,  ^  ,  2 )  =  g  {x,l,z)/h  {x  ,  z, 
=  F(x,i).  If  t'?)  or  r.(z)  is  a  rational  function  of  2,  then  the 
more  useful  expressions  of  g"*”  (x,?  ,z)  and  g  (x,^,z)  have  been 
obtained  in  1211  (equations  43,  44,  50  and  51) • 

(3J.6  )  can  now  be  re-arranged  as 

(5.13)  g"*"  (x ,  C  ,  z)  ^(x,5,z)  -  (x  ,  ?  ,  z )  exp  (-zw)  } 

+  xj^f  g  (x,  ^  ,z)  exp  (-zD)  n  (' -z)  ^(x,-;,z)} 

=  x™''’^T*{9  (x,C,2)exp(-zw)  } 

-  XT*''  ?  (x,- ,2)exp(-zD)Q('-z)Qj^^^(x,' ,z)  } 
xg"(x,^,z)Cjjj^^(x,5  ,-z) 

for  0<|x|sl,  Re(C) ^Re(z) 2O ,  m>0,  w>0.  As  the  left  hand  side  of 
this  equation  belongs  to  and  its  right  hand  side  belongs  to 
Ng/  applications  of  Liouville's  Theorem  and  analytic  continuation 
yields 

(5.14)  g''’(x,C,z)W^  ^(x,C,z)  -  x‘'”’^^T2^5  (x  ,  C  ,  z)  exp  (-zw)  } 

xTg^  9  (x,C,z)exp(-zD)n(C-z)Qjj^  ^(x,C,z)  }  =R(x,C) 

for  0<|x|<l,  Re(4)^0,  Re(z)>0  ^  m>0,  w>0;  and 

(5.15)  x"''’'^Tjt  g  (x,?,z)e:q>(-2w) }  g"(x,' ,z)exp(-2D)n(F-z)Q^  ^(x,5,2) } 

-  xg“(x,C,z)Cj^^^(x,C,-z)  ®  R(x,?) 

for  0<|x|sl,  Re(5)20,  Re(z)^0,  mjO,  wjO. 


_  T  a  _ 


Also,  observing  onao  iron  ;3.12', 


W  (x,  1 , 2)  exo  (zD'  1  (z)  -  Q  i:<,L,z) 
n ,  w  -  n ,  w  ' 

i  -—1  k  I 

=  E|ik=in+1  (^in(Wj^+Sj^-D,0)  )  )  ,|nQ=in,WQ=w  j  < 

we  re-arrange  (3.16)  differently  and  then  apply  Liouville's 
theoren  and  analytic  continuation  to  obtain 


N_ 


(5.16) 


Q  (x, ' , z) 
m,w  ^ 

g"*"  (x,  ^  ,  z) 


jn+i  [  exp(2  (D-w)  )'M2)  ] 


-  X  j  ^ 

'  z ! 


g  (x, 1 , z) 


f  I 

(x,?,-z)exp(2D)'i' (2) 

*  V,  ; 


S  ( X  ,  2  ) 


g  (x,C/Z) 


i 


for  0<|x|<l,  Re(C)>0,  Re(z)>0,  m>0 ,  w>0 ;  and 

/^  /s 

W_  ,  (x,  5  ,  z)  exp  (  zD) (z)  -  Q  (x,C,z) 


(5.17)  - 


n,  w 


m,  w 


g“ (x, C  »  z) 


f  exp  (z  (D-w)  )  4*  (z) 

■  - 

C  (x,C,-z)exp(zD)  ■f' (z)  ] 

m ,  w  -  ( 

g'*’  (x,  z) 

g'''(x,C,2)  1 

J 

1 

m+1  . 

+  X 


for  0<lx|.^l,  Re(O>0,  Re{z)<0,  m>0,  w^O , 


S  ( X  /  ^ ) 


The  expressions  of  R(x,C)  and  S(x,5),  which  are  dependent 

on  the  type  of  busy  cycle  we  are  interested  in,  will  enable  us 

to  find  the  expressions  of  0  (x,C/2)  and  C  (x,5,-2)  as  in 

in ,  w  m ,  w 

the  following  theorems: 

THEOREM  7  :  C®  (x,6,-z)  and  of  ,  (x,^,z)  satisfy  the  following 

— — —  m,w  in,w 


simultaneous  integral  equations; 


■  -2^ , 


(5.19) 


X  ,  ;  ,  2  ) 


(  O'  X  1  j  Rs  ! ; !  >.j  ,  Rs  '.  z  ^  3  ,  ,  v.’^O 

m-1 


-  X 


+  xc“  ./'X,  -: ,  0' 
m  /  ^ 


:  X ,  1 , 2 . 


^  m+l, 
X  T. 


exp ( 2 (D-w) ) M  2) 
g* (X, ^ , 2) 


-  XT  ^ 

'  z 


g  { X ,  -: ,  3 ) 

^  s 

,(X  ,  ,  -2)  exp  (  zD  '  •  z  . 
:  m ,  _ _ 


c  ( X  ,  C  ,  2  ) 


(0<;x4l,  Re^-:)>0,  Re;2;>0,  zi>0,  w>0.' 


PROOF:  We  have  p{i^  =  0 '■  =  0  for  all  v>l .  Hence  li.T-  [y.,1. 

Thus  if  we  let  2-*-=  in  (5.15),  we  shall  obtain 

(5.20)  R  {x,^]  -  X®  ^li.Ti  (x,  C  /  z)  exp  (-zw)  } 

-  xli.Tij,^^^  T*{  g  (x,' ,2)exp(-2D)  :.  (;-z)Q^^,_^(x,^  ,2)  } 

Upon  substituting  this  back  into  (5.15),  on  behalf  of  (5.3), 

we  obtain  (5.18).  Also, we  have  lim^^^Qg'*’ (x,  ■  ,  z)  ^0  and  from  (D.43), 

lim  (x,C,z)  =  0.  Thus  if  we  let  z-0  in  (5.16),  remembering 

that  lim_  _  T  =  lim  ,,}(z),  we  obtain 

z-^o  z  Z*(J 


(5.21)  S(x,5) 


-  X 


m+1 


xC  (x , 
m,w 


0) 


g'*'(x,5,0) 


Upon  substituting  this  back  into  (5.16) 
we  obtain  (5.19).  Q 


simultaneous  integral  equations; 


(5.22)  g' (.x  ,  t ,  z  )  (x  ,  , -z)  =  x"'v*ig  z)  e:<p  i-zw)  } 

in  f  v\  z 

-  v*^g  (X,  ^  ,  Z)  exp  (-ZD)  .-  (  i-Z  )  (x  ,  ,  7  )  ; 

(  0<ix!.<l,  ReC.^^O,  Re(z)<^0,  m>0,  w>0) 

(5.23)  q  (x,i,z)CQQ(x,-;,-z)  =g  {x,^,z)  -  1  im^  _  ( x  ,  ■  ,  z  ) 

-  V*(g  (x, -;,z)exp(-zD)  .1  (;-z)Qq  q(x,{,z)} 
(0<|xl<l,  Re(S)>0,  Re{z)<0) 


(5.24) 


g"*^  (x,  ?  ,  z) 


-  x”-^^  a  xc;;^^(x,-  ,0) 
g'*’  (x,  5 , 0) 


exp  (z  (D-w)  )  f  (z) 

+  x  T  •  - - - 

z  +  <  r  \ 

g  x,C/Z) 


Cm^„(x,C  ,-z)exp(zD)  Hz) 
g'*’  (x,  ?  ,  z) 


(0<jxl<l,  Re(C)>0,  Re(z)>0,  m^O,  w>0) 


PROOF !  We  have  p{w^=0}  =  0  for  0<mQ<l<^nj^.  Hence 

/N 

~  0  where  6^  ^  is  the  Kronec)<er  delta. 

Thus,  for  m>0,  w>0,  letting  z-*-®  in  (5.14),  we  obtain 

(5.25)  r'^Cx,?)  =  -x"’^^lim  j  {g  (x ,  C  ,  z )  exp  (-zw)  } 

2  *♦•00  *2 

+xliin  T,f  g~(x,?,z)exp(-zD)0(^-z)Q!'  (x,C,z)  } 

Z  Z  Til  f  W 

Upon  substituting  this  baclt  into  (5.15),  on  behalf  of  (5.5),  we 
obtain  (5.22).  When  mQ=0  and  Wq=0,  we  first  modify  (5.14)  to  have 
the  term  q(x,5,z)-:^  included.  (5.23)  can  now  be  derived  by  the  same 

metnod  as  that  for  (5.22).  For  (5,24),  the  proof  is  similar  to 
that  for  (5.19) .  Q 


\ 


Rer.arks  : 

,?,  .6]  When  ,  [5.  IS)  becomes 

(5.2oj  C^^^(x,C,-z)  =  x‘’^(j*{g  (X,  ;  ,  z)  exp  {-zw)  ;■  /  g~(x,;,2) 

which  is  equation  (198)  in  [24] .  (3.16)  now  can  be  written  as 

(  =  •27)  '''ni,w^^' (Jz'--  ^  '  2)  exp  (-zw:  ;/g  ■  (x  ,  i  ,  2)  . 

This  equation,  together  with  (3.1),  (3.2),  (5.6),  (5.26)  yields 

(5.23)  ^(x,^,z)  =  {g  (x , '  ,  z)  ex?  ( -zw) //g"^  (x ,  ' ,  z )  , 

/ 

a  well-known  result  due  to  Pollaczek  [14] ,  Kingman  [10]  and  Takacs 

[21]  . 
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PROOF :  Ir  -  x-exp(--t)  for  t>0,  then  for  all  v>l, 

is  independent  of  both  and  ^  and  E '  exp  ( -zi,^)  }  =  \/ {',  +  z) 

for  Re(z)^0.  Thus  we  obtain  (6.1).  (3.16)  can  now  be  written  as 

(6-5)  W^^^{x,l,2)  = 

m+l  ,  y  ^  .  e  ^  ^ 

X  (X  +  --Z)expi-2W)  -  ,\XP=  ^(x,l,0)  -  \xexp(-zD)Q®  (x,l,z) 

— - - - _ _ _ m ,  w 

A  +  '  -  2  -  A  X  7  (  2  ) 

for  0<|xj^l,  Re(4)>0,  Re(2)i.O,  m^O,  w^O.  Now,  since  W®  (x,',z) 

in ,  w  ^ 

is  analytic  in  the  domain  Re(2)^0,  letting  z=e  as  defined  in  (D.44) 
yields 


(6.6)  P„,„(x,?,0)  - 


-9D"s 


for  0<|x|<l,  Re(')^0,  mjO,  w>0.  Also,  from  (5.19),  we  can  write 


(6.7) 


g'^(x,C/z)  g'^(x,C,0) 


m+i  exp(z(D-w) )7 (z) 

+  X  T  ‘  - —  > 

2  +  ,  r  ^  I 

g  (X,CrZ) 


exp  (zD)  'i'  (z) 


+  AxP^  ^(x,c,o)  — ^ - -  T 

g  (x,C,0)  (A+C)  ^  g''’(x,C,z)  (A+5-Z)  j 

for  0<|x|.$l,  Re(^)j.O,  Re(2)^0,  m>0,  w>0.  As  g'*’(x,4,z)  ta)ces  the 
form  of  (6.4),  we  let  z-^9  in  (6.7)  and  then  eliminate  Q®  ^(x,^,6) 

from  the  resulting  equation  and  (6.6)  to  prove  (6.2).  Q 
Remarlc ; 

(R.6)  When  ,  (6.2)  and  (6.6)  will  become 

a  well-known  result  for  the  residual  busy  period  of  the  M/G/1  queue  ; 


(0<|kUi,  Re(C)>0,  m>0,  w>0) 
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GI/G/1  queue.  Limited  virtual  waiting  times.  Finite  dam. 
Regenerative  processes.  .Actual  waiting  times.  Lost  service 
times.  Virtual  waiting  times.  Busy  cycle.  Transient  results 
Stationary  results. 
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This  paper  studies  the  GI/G/1  queueing  system  in  which  no 
customer  can  stay  longer  than  a  fixed  interval  D.  This  is  also 
a  model  for  the  dam  with  finite  capacity,  instantaneous  water 
supply  and  constant  release  rule.  Using  analytical  method 
together  with  the  property  that  the  queueing  process  "starts 
anew"  probabilistically  whenever  an  arriving  customer  initiates 
a  busy  period,  we  obtain  various  transient  and  stationary 
results  for  the  system. 
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